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Abstract. The current status of NLTE model atmosphere calculations of A type stars is re- 
viewed. During the last decade the research has concentrated on solving the restricted NLTE line 
formation problem for trace elements assuming LTE model atmospheres. There is a general lack 
of calculated NLTE line blanketed model atmospheres for A type stars, despite the availability 
of powerful methods and computer codes that are able to solve this task. Some directions for 
future model atmosphere research are suggested. 
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1. Introduction 

Compared to the atmospheres of hotter O and B stars that have stellar winds and to 
atmospheres of cooler F and G stars that have convective atmospheres, the atmospheres 
of A-type stars are relatively quiet, which enables the presence of various interesting 
phenomena. This calmness of atmospheres leads to the development of a number of 
different chemical peculiarities and magnetic field structures on long time scales. 

There exists a widely used grid of line blanketed LTE model atmospheres fKurucz ll993l 
Castelli & Kurucz f2(JU3fl . which is also extensively used for A stars. Besides this grid, 
several attempts to remove the inconsistent assumption of LTE have been performed. 
The very first NLTE calculation of a cool B star, which can be somehow understood as a 
first step toward A star modelling, was performed by Auer & Mihalas p97l)fl . The domain 
of A star temperatures was first reached by Kudritzki H1973(l and Frandsen H1974|l . A 
grid of NLTE model atmospheres of stars with 10000 K < Tos < 15000 K (thus also 
including a part of the A star domain) was calculated by Borsenberger & Gros H1978|l . 
A great stride was made by Hubeny l|1980l I1981|l . who studied in detail the effect of 
the La line on the NLTE model atmospheres of A stars and applied his code to star 
Vega. A detailed analysis of the physics which is contained in releasing the assumption 
of LTE was presented by Hubeny H1986|l . First NLTE line blanketed model atmospheres 
of A stars were calculated using a method of superlevels and superlines by Hubeny & 
Lanz p993|l . Unfortunately, although there was great progress in calculating NLTE model 
atmospheres for hot stars and for white dwarfs, calculations of NLTE model atmospheres 
of A stars are very rare. The only attempt of a NLTE model atmosphere of an A-type 
star was the model calculated by Hauschildt et al. (1999j for parameters corresponding 
to Vega. The authors announced a grid of NLTE models to be published, but till now 
nothing appeared. 

The inherent reason which stands tacitly behind the lack of calculated NLTE A star 
model atmospheres is in the different scaling properties of the atmospheres of this stellar 
type. The atmosphere of A stars have very extended line formation regions meaning that 
different lines form at very different depths (cf. Fig. 1 in Lanz & Hubeny, 1993). Therefore 
they are very difficult to describe by a simple standard depth discretization as compared 
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to O stars and white dwarfs, and hence a much more elaborate depth scale has to be 
used. Our own experience shows that there are severe convergence problems caused by 
improper depth scaling and only a tremendous number of depth points helped. 

2. Basic equations for NLTE model atmospheres 

A solution of a model atmosphere problem means taking the basic input values of 
stellar luminosity L, stellar radius R, and stellar mass AI (or, equivalently for the usual 
plane-parallel case effective temperature Tcff and gravitational acceleration at the stellar 
surface g), adding physical laws that are important in the stellar atmosphere (a principal 
ingredient), and calculating the space distribution of various physical quantities (tem- 
perature T{r), electron density ne(r), population numbers ni(f), density p(r), velocity 
field vif), etc.). This procedure is often being performed in substeps. The first step is the 
determination of T, n^, p, and rii for the most important opacity sources. The second 
(optional) step is the determination of rii for the minor abundant species, often referred 
to as trace elements. The final step in this process is the calculation of the theoretical 
emergent radiation which may be then compared to the observed one. 

A standard model atmosphere here is a static one-dimensional (plane-parallel or spher- 
ically symmetric) atmosphere in hydrostatic, radiative, and statistical equilibrium. To 
calculate such a model atmosphere we need to solve the radiative transfer equation (to 
determine the radiation field /), the hydrostatic equilibrium equation (to determine p), 
the energy equilibrium equation (to determine T), and the statistical equilibrium equa- 
tions (to determine rii). 

2.1. Radiative transfer equation 

The radiative transfer equation is the principal equation for model atmosphere calcula- 
tions. Regardless of the geometrical approximation which is adopted, the basic solution 
of the radiative transfer equation is performed along a ray, which reads 



Usually the simplest possibility, i.e. the plane-parallel atmosphere, is being used for the 
modelling of stellar atmospheres, basically due to its relative simplicity. 

2.1.1. Formal solution 

By formal solution of Equation H2.1|l we mean the solution for a given opacity and 
emissivity. Such a solution may be performed using either differential or integral methods. 
The latter are especially useful for the multidimensional case. The inconsistency of the 
formal solution is that opacity Xi^ and emissivity rj^ are not given, but in the process 
of solving the model atmosphere problem they depend on temperature, density, and 
radiation field, which are to be determined when solving the radiative transfer equation. 
However, the formal solution remains at the heart of each model atmosphere code (see 
Auer lMT^ . 

2.1.2. Approximate solution using ALO 

The process of the formal solution may be expressed using the A-operator 



d/.(r) 
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(2.1) 




(2.2) 



This expression may be used in an iterative process, where the source function Si, is 
determined by solving the other constraint equations of hydrostatic, energy, and statis- 
tical equilibria. As has been shown, such a process, however effective for the optically 
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thin cases of molecular clouds (cf. Dickel & Auer ll9iJl|l . fails to produce a convergent 
solution for the optically thick cases of stellar atmospheres (see Auer ll98^ . This problem 
has been overcome using the Newton Raphson method for model stellar atmospheres by 
Auer & Mihalas H1969|l . As an efficient alternative, an approximate A-operator A* may 
be used, 

J, = A*5, + (A, - a;) (2.3) 

The operator A* is constructed in such a way that it contains the basic physics of the 
problem and is simply calculable. The most efficient method is to calculate this oper- 
ator consistently with the numerical method used for the formal solution, like it has 
been described by Rybicki & Hummer H1991|l or Puis H1991|l . Then the radiation field is 
determined by an iterative process, 

4") = a: [s, (4"))] + (A, - A^) [S, (2.4) 

where (n) means the current iteration step and {n — 1) the previous one. 

2.2. Energy equilibrium 

The energy equilibrium equations determine the temperature structure. They describe 
the basic energy balance in the atmosphere. For the case of static atmospheres, the 
equation of radiative equilibrium, 

/•oo 

{k^Jjj - Tjv) dv = I {Jv - Sv) dv = 0, (2.5) 

which simply describes the radiative energy balance, is usually being used. However, this 
equation does not guarantee the radiative flux conservation, which results in an incorrect 
temperature structure at large optical depths. Therefore, another form, which comes 
from V • = 0, is to be used at large optical depths. Better numerical properties of the 
scheme may be achieved if a linear combination of these equations is considered, as was 
first done by Hubeny & Lanz (jl995|l . Another improvement of convergence properties, 
especially useful at low continuum optical depths where strong lines are still present, was 
introduced by Kubat et al. where a thermal balance of electrons instead of the 

radiative equilibrium equation is used. 

Then three different equations may be used in the stellar atmosphere to take into 
account energy equilibrium. In the deepest layers, the differential form of radiative equi- 
librium is used. In the middle parts, where continuum radiation is formed, the integral 
form of radiative equilibrium works well. In the outer parts where optically thick lines 
coexist with an optically thin continuum (which is the case of A star atmospheres), the 
thermal balance of electrons works best. A more detailed discussion with the numerical 
properties of these methods may be fomid in Kubat H2003all2003b|l . 



2.3. Equations of statistical equilibrium 

The equations of statistical equilibrium are the key equations both for the NLTE model 
atmosphere and line formation problems. The set of equations for statistical equilibrium 
for the static case may be written as 

Y^^Ra +Ca)~Y. "'(-^'^ + Cu) ^ (2.6) 

I I 

Radiative rates Ru are those responsible for the NLTE effects, since they introduce non- 
local interaction and they alter the level population irrespective of local equilibrium con- 
ditions. They move the gas outside thermodynamic equilibrium. On the other hand, for 
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an equilibrium (Maxwellian) electron velocity distribution (which is the common case), 
collisional rates reintroduce thermodynamic equilibrium. The balance between coUisional 
and radiative rates determines the applicability of the LTE approximation. If collisions 
dominate, then the LTE approximation is feasible. If radiative transitions dominate, 
which is the case in A star atmospheres, then equations of statistical equilibrium need 
to be solved if we want to obtain the correct population numbers. 

Since the system of rate equations is linearly dependent, we have to close it with 
some other condition. For a reference (usually, but not necessarily hydrogen) atom ei- 
ther particle conservation equation J^k — N — Ue or charge conservation equation 
YlkTljlj^jk = ne are used. For other atoms, the abundance equation Nk = YkNr is 
usually used. Detailed forms of the equations may be found, e.g., in Kubat 1)1997(1 . 

2.4. Solution of the system of equations 

The whole system of equations is usually solved using a Newton-Raphson iteration 
scheme, which in modeling stellar atmospheres used to be referred to as linearization 
( Auer & Mihalas 11969(1 . The radiation field is linearized as well, or it can be included 
using approximate lambda operators (see Eq. 12.4(1 . as has been first done by Werner 
(WM). Excellent reviews of ALI methods of model atmosphere solutions are provided by 
Hubeny (,2003;1 and Werner et al. I(2()()8(l . In addition, an implicit linearization of 6- factors 
saves additional computer time and memory f Anderson 11987(1 . The whole process may 
be accelerated using either Kantorovich (see Hubeny & Lanz 1193^ or Ng acceleration 
(see AuerHnmi. 

2.5. Line blanketing 

Line blanketing is caused by enormous amounts of spectral lines in the UV region (~ lO''), 
mostly of iron and nickel. This enormous amount of lines causes radiation to be absorbed 
in the UV and reemitted in the visible region, thus changing the emergent flux and 
atmospheric structure significantly. 

There are two main approaches to line blanketing under the simplifying assumption of 
LTE, the ODF (Opacity Distribution Function) and OS (Opacity Sampling). 

If we do not use the assumption of LTE, we have to solve the equations of statistical 
equilibrium to obtain the correct population numbers for all levels. This task is relatively 
simple for simple atoms (like He), but becomes difficult for atoms with a complicated 
level structure, like iron. In order to cope with the complexity of these atoms, Anderson 
1(1989(1 introduced the concept of superlevels and superlines. A superlevel is a level which 
is created by grouping several individual levels. In order to achieve a real simplification, 
it is useful to group the levels in such a manner that the relative population distribution 
inside the superlevel obeys the LTE distribution. Thus, this grouping is done according 
to the energy, Ei, of individual levels (Anderson'1989, Dreizler & Werner 1993"). A more 
sophisticated grouping according to Ei and parity was done by Hubeny & Lanz 1(1995(1 . 
Dreizler & Werner 1(1993(1 used opacity sampling in their calculations, whereas Hubeny 
& Lanz 1(1995(1 used the NLTE opacity distribution function. 

2.6. LTE versus NLTE 
Since the beginning of the NLTE calculations there is a battle between 'LTE people' 
and 'NLTE people', which of the approaches is better. If NLTE model atmosphere is 
calculated as easily as the LTE one, there would be no discussion and everybody would 
use the NLTE one. However, this is not the case and it is much more difficult to calculate 
a NLTE model than the LTE model. Two conflicting aims enter the scene. The first one 
is to analyse as many stars as possible, which can be hardly done using NLTE models. 
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The second one is to study stars as accurately as possible, which can be hardly done 
using LTE models. 

For the case of LTE it is relatively easy to handle line blanketing, since one neglects 
the effect of radiation on atomic population numbers, which saves computing time enor- 
mously. On the other hand for a more general case of NLTE, including line blanketing 
is computationally expensive, albeit nowadays it is possible to handle NLTE calcula- 
tions using sophisticated numerical methods and contemporary computers. However, for 
extended atmospheres the reliability of LTE decreases and one is forced to use NLTE. 

To avoid repeated calculations of complicated model atmospheres one may take the 
advantage of precomputed grids. However, one has to keep in mind that using such grids, 
which were calculated under certain physical assumptions, always limits the results by 
these assumptions. 



3. NLTE line formation calculations 

Calculations of full NLTE model atmospheres of A type stars is a difficult task. Due 
to a large span of line formation regions down to optical depths in the continuum of 
about 10^^, the numerical procedure becomes extremely unstable and calculations need 
special care. That is why people started to solve easier task of NLTE line formation for a 
given model atmosphere. A necessary condition for such calculation to be reasonable, is 
negligible inffuence of the ion on the ionization balance and, consequently, on the global 
structure of the atmosphere. Such elements are being referred to as trace elements. 

Contemporary NLTE analysis of A stars are mainly concerned with solving the statisti- 
cal equilibrium equations for trace elements for a given (mostly LTE) model atmosphere. 
These kind of calculations have already been reviewed by Hubeny (1986^ and Hubeny 
& Lanz (1993), so only new calculations that appeared after the last review are listed 
in Tabled We did our best to mention all calculations. In the case we unintentionally 
omitted some, we apologize for it. 

There are several codes available for this purpose. They use slightly different numerical 
techniques, basic differences are in the treatment of the radiation field. Both the accel- 
erated lambda iteration and complete linearization methods are used. A comparison of 
results from different codes for Vega (and Sun) was presented by Kamp et al. H2(J(J3|I . 
who found large differences in the results different codes. Therefore we added to Tabled 
an indication which code was used. 

4. Beyond standard model atmospheres 

4.1. Nonthermal collisional rates 

In standard NLTE calculations only the radiation field and level populations are allowed 
to deviate from their equilibrium values. The velocity distribution of individual species 
is still assumed to be in equilibrium. Collisional rates are then calculated as an average 
over the Maxwellian electron velocity distribution. 



Thus in equilibrium there holds a relation between collisional rates up and down (exci- 
tation/deexcitation or ionization/recombination). 




(4.1) 



(4.2) 
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Table 1. List of NLTE calculations for a given LTE model atmosphere for A stars 
Ion reference code comment on included levels and transitions 



He I Takeda flTOIl 3 

Li I Mashonkina et al. 1200211 4 

Lil/ll Shavrina et al. I2UU11 

C I Taked a (,199 2b)l 3 

Venn lP??K)l 2 

Rentzsch-Holm l|1996b|l 2 

Paunzen et al. (1999) 2 

Cl/ll Przybill a et al. L2001bJ 1 

Nl Takeda H1992b|l 3 

Sadakane et al. 11993(1 3 

Takeda & Takada-Hidai (tlM^ 3 

Takad a-Hid ai & Takeda (IB^ 3 

Venn iP??^ 2 

Lemke & Venn fT5M^ 2 

Rentzsch-Holm (1996a) 2 

Nl/ll Przybilla & Butler (2001J 1 

Ol Takeda (1992a, 1997J 3 

Takeda & Takada-Hidai llB^ 3 

Takeda et al. ( 199^ 3 

Paunzen et al. ( 19^ 2 

Przybilla et al. ^2W)^ 1 

Nal Takeda & Takada-Hidai (pMIji 3 



Mashonkina al. H2000|l 4 

Mgl Shimanskaya et al. 1 2000 t 4 

Idiart & Theve nin ifgOOC^ 5 

Mgl/ll Przybilla et al. 12001all 1 

S I Takeda & Takada-Hidai (IT995)l 3 

Takada-Hidai & Takeda (TPB:) 3 

Kl Ivanova & Shimanskii (;2(30i) ) 4 

Cai Idiart fc Thev enin (:'2000.) 5 

Till Becker ifUM)! 1 

Fell Becker llT^^ 1 

Fel/ll Rentzsch-Holm (fl996b| 2 

Thevenin & Idia rt il999i 5 

Srll Belyakova et al. H1999|l 4 

Ndll/lll Mashonkina et al. 1200411 1 



calculations for Deneb; 88 levels 
20 levels 

129 levels, 2351 radiative transitions 
83 levels 

83 levels, 63 transitions 

83 levels, 63 transitions 

C I levels with n ^ 9, C ll levels with n ^ 4 

119 levels, 2119 radiative tra nsitions 

used model atom of Takeda (11992b I 

used model atom of Takeda ( 19^^ I 

used model atom of Takeda (Il992b I 

93 levels, 189 transitions 

93 levels, 189 radiative transitions 

96 levels, 82 transitions 

N I levels with n ^ 7, N ll levels with n ^ 6 

86 levels, 294 radiative transitions 

86 levels, 294 radiative transitions 

calculations from Takeda l|1997|l for CP stars 

15 levels 

all levels below excitation energy 
calculations for A supergiants, 92 levels, 178 
radiative transitions 

solution for T^a = 4000 - 12 000 K, 21 levels 

solution for T^a = 4500 - 12000 K, 50 levels 

104 levels, 980 radiative transitions 

all levels with n ^ 10 

56 levels, 173 radiative transitions 

56 levels, 173 radiative transitions 

solution for Teff = 4000 - 10 000 K, 36 levels 

84 levels, 483 transitions 

complete model atom, using superlevels 
complete model atom, using superlevels 
79+20 levels, 524-23 lines 
2564-190 levels, 2117+3443 lines 
solution for T^a = 4000 - 12 000 K, 41 levels 
247+69 levels 



Computer codes used for calculations: 1 - DETAIL (^Ciddines I1981|l : 2 - Kiel (Steenbock & 
Holweger[THHU; 3 - Takeda HWHi : 4 - N0NLTE3 fSakhibuUin IBH^ : 5 - MULTI (Carlsson 



If the velocity distribution is not Maxwellian (e.g. in electron beams formed in the Sun 
during flares - cf. Kasparova & Heinzel J0()2J, then the equation (|4.2|> . which expresses 
the equiUbrium condition, is not valid. In such a case the collisions may cause level 
population numbers to differ from the LTE values, which has an observable effect on the 
line profiles. The nonthermal coUisional term may be considered as another source of 
NLTE effects. Such effects are present in the Sun, where the electron beams are formed 
after magnetic reconnection. There is also a possibility that they may be present in 
magnetic A stars, where releasing magnetic energy in eruptive events like flares may be 
expected as well. 
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4.2. Magnetic fields 

Inclusion of magnetic fields into model atmosphere calculations have been done only 
occasionally. An important attempt in calculating a model atmosphere with a magnetic 
field was done by Carpenter H1985|l . who found changes in the net gravity and pressure 
distribution due to the magnetic field. The most consistent model so far was recently 
developed by Valyavin et al. H2(JU4|I . who assumed LTE and included the Lorentz force 
into the equation of hydrostatic equilibrium. A historical overview of magnetic model 
atmosphere calculations is also presented there. 

4.3. Diffusion 

An important aspect of A star atmospheres is the fact that atmospheres of a good fraction 
of stars of this stellar spectral type are extremely quiet. Such quiet atmospheres without 
significant global motions like global convection and stellar winds enable long time scale 
processes of diffusion and gravitational settling to take place. Diffusion occurs not only 
in the presence of magnetic field (note the effect of enhancing radiative acceleration in 
polarized radiation field and the effect of ambipolar diffusion, Babel & Michaud I1991al 
I1991b|l . but also due to different sensitivities of various atoms and ions to incoming 
radiation. Radiative diffusion is able to explain various chemical peculiarities in CP 
stars and is also responsible for the isotopic shift effect (Aret & Saoar f^02|l . However, 
consistent NLTE model atmospheres that take into account radiative diffusion processes 
are still missing for A type stars. NLTE model atmospheres with diffusion were calculated 
by Dreizler & Shuh H2l)U3|l for white dwarfs, where diffusion processes are also important. 
The physical process of diffusion is reviewed in detail elsewhere in these proceedings 
(Michaud). 

4.4. Extended atmospheres 

Plane-parallel model atmospheres are also being used for NLTE modeling of atmospheres 
of A-type supergiants (e.g., Kudritzki I1973|l . However, they do not describe the limb 
darkening correctly. A-type supergiants have also stellar winds and show P Cyg line 
profiles, therefore it is impossible to describe them using plane parallel atmospheres. 
Also, wind models of A stars are not calculated too often. Multicomponent hydrodynamic 
radiatively driven wind models for A stars were calculated by Babel H1995|l . Another 
radiation driven wind model was calculated by Achmad et al. H1997|l . but they did not 
take into account non-LTE effects. Generally, for extended atmospheres the NLTE effects 
are stronger. An attempt to calculate spherically symmetric NLTE model atmospheres 
with a stellar wind of Deneb was done by Aufdenberg et al. I|2()()2|l . However, as the 
authors note, they were not able to fit the observed Ha line profile. Work on developing 
consistent NLTE wind models of A supergiants is currently in progress (Krticka & Kubat, 
these proceedings). 

4.5. Limh darkening 

An important property of emergent radiation from stellar atmospheres is the angle de- 
pendence of the specific intensity /(0), usually called limb darkening. Approximate limb 
darkening laws are often used (see, e.g., Allen. [TI)63I or Grav. IT^THIl . These approximate 
laws don't describe the angular dependence of the specific intensity very well even in 
the case of thin stellar atmospheres. In addition, limb darkening is strongly frequency 
dependent (Hadrava & Kubat 2003), and in the center of a line, even limb brightening 
instead of darkening may appear, as can be seen in Figure^ This effect becomes very 
strong especially for extended stellar atmospheres. 

For a correct description of limb darkening it is necessary to use model atmospheres 
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Figure 1. The comparison of limb darlcening for the thin atmosphere (left panel) and the 
extended atmosphere (right panel). These results were obtained from the Korcakova & Kubat 
code lIMiHl . 

with a better geometry than plane-parallel. This has been done by Claret & Hauschildt 
(,20021), who, using a spherically symmetric code, calculated limb darkening for a grid of 
model atmospheres from A to G spectral types. 

Knowledge of an accurate limb darkening law is necessary for detecting stellar spots 
(see Kjurkchieva, il989p . It is also very important for the analysis of interferometric 
measurements, as well as for the correct treatment of stellar rotation. Luckily, eclipses 
from binaries (e.g., Twigg [T979|l . and the gravitational microlensing effect (Heyrovsky 
I200H|) measurements of limb darkening are possible for at least some of the distant stars. 

4.6. Rotation 

Stellar rotation is usually being accounted by convolving the rotation profile with the 
static line profile (see Grav ll976p . This technique is inaccurate not only due to using only 
approximate limb darkening laws, but also because the dependence of limb darkening on 
frequency as well as gravity darkening are neglected. A more accurate method is to 
integrate the static intensity over the rotating disc of the star, where gravity darkening 
can also be included. A "simple" way of including gravity darkening in this integration is 
described by Collins 11963p . This method was used by Gulliver et al. ()1994ll to determine 
the rotation of Vega. 

Stellar rotation is not like a rigid body rotation, as it is usually assumed, but it is 
probably differential. Until now, there are only few measurements of differential rotation. 
Unfortunately, in many works, the authors use the simplest form of the limb darkening 
law and the obtained results can be used only with reservations. 

For an exact description of stellar rotation the plane-parallel model atmosphere is 
insufficient. It is necessary to use a hydrodynamic multidimensional model atmosphere, 
at least an axially symmetric one together with the correct description of the radiation 
field. Recently, we started to work on such a model (Korcakova et al, these proceedings). 



5. Conclusions 

Full NLTE model atmospheres of A stars are difficult to calculate, since the lines of 
different ions form at very different depths, and therefore there is a need for a huge 
amount of depth points to resolve all line formation regions sufficiently. Only LTE model 
atmospheres together with NLTE calculations for individual ions have been calculated 
in significant amounts. This underlines the need for more intensive calculations of NLTE 
model atmospheres of A stars. Adding consistently new different physical processes to 
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NLTE model atmospheres calculations, such as diffusion (see Michaud, these proceed- 
ings), magnetic fields (see Moss, these proceedings), and more accurate treatment of 
convection zones (see Kupka, these proceedings), is necessary. This is the real challeng- 
ing task the near future. 
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